Quantum chaotic resonances from short periodic orbits 
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We present an approach to calculating the quantum resonances and resonance wave functions of 
chaotic scattering systems, based on the construction of states localized on classical periodic orbits 
and adapted to the dynamics. Typically only a few of such states are necessary for constructing a 
resonance. Using only short orbits (with periods up to the Ehrenfest time), we obtain approximations 
to the longest living states, avoiding computation of the background of short living states. This 
makes our approach considerably more efficient than previous ones. The number of long lived states 
produced within our formulation is in agreement with the fractal Weyl law conjectured recently in 
this setting. We confirm the accuracy of the approximations using the open quantum baker map as 
an example. 
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Quantum scattering in chaotic systems is a very active 
field in the area of quantum chaos, with current experi- 
mental realizations in ballistic quantum dots pj], micro- 
lasers [1] and microwave cavities among others [1]. 
The fractal Weyl law f|, which relates the counting of 
resonances in the complex plane to the dimension of the 
trapped set of the corresponding classical dynamics, has 
attracted considerable attention 0, 0, H, [1] because reso- 
nances (or Gamow states) are central to the description 
of many aspects of wave scattering. 

The decaying eigenstates associated with these quan- 
tum chaotic resonances are far from being fully under- 
stood. They have recently been shown [1, [l^] to display 
fractal structures in phase space when the resonance is 
long-lived, that is when the decay rate T/h remains fi- 
nite as ?i — !■ 0, and to be localized when the resonance is 
short-lived, that is when ft/T — s- 0. However, the semi- 
classical limit is much richer for scattering systems than 
for closed ones, for which the quantum ergodicity the- 
orem ^11] states that almost all states become uniform. 
Owing to the existence of different decay rates, nothing 
of this kind is available for scattering systems and we are 
still far from a complete description. 

Our purpose here is to establish an approach to res- 
onances and resonance wave functions based on short 
classical periodic orbits. The idea is to use the prolif- 
eration of periodic orbits in the phase space of chaotic 
systems to build an approximate basis of functions for the 
quantum Hilbert space. These functions are constructed 
in such a way as to contain dynamical information up 
to Ehrenfest time. This formulation has several virtues. 
First, the fractal Weyl law emerges very naturally from 
the theory and is seen to have a direct connection with 
periodic orbits. Second, we have an approximation to the 
quantum propagator that provides the long-lived states 
(which are usually dominant) without having to calcu- 
late short-lived states, therefore very significantly reduc- 
ing the dimension of the matrices involved in the theory. 
Specifically, we achieve a power saving in the matrix di- 



mension. Third, it turns out that usually only a few of 
our states are required to produce a quantum resonance, 
providing a way to quantitatively analyze scarring effects 
(anomalous localization of chaotic eigenstates around pe- 
riodic orbits fl^). Finally, it opens a new and promising 
avenue for semiclassical approaches to resonance wave 
functions, which have so far been elusive. 

A corresponding theory exists for closed systems in the 
form of scar functions [isl . [T^ . which have proved effi- 
cient in providing semiclassical approximations for quan- 
tum spectra and eigenstates of billiards and quan- 
tum maps In the open systems considered here, 
the efficiency gain is considerably greater. An alterna- 
tive periodic orbit approach to resonances already exists 
in terms of the semiclassical zeta function [l^. However, 
the orbits used are in general much longer than the ones 
considered here and the final result is an approximation 
to the spectral determinant that does not provide the 
wave functions. 

For simplicity, we restrict ourselves to quantum maps, 
in which time evolution is discrete, the quantum Hilbert 
space has finite dimension N = l/(27r/i) and the classical 
phase space is a torus. Open maps are defined by iden- 
tifying a region of phase space -usually a strip of width 
M /N parallel to one of the axis- with a 'hole', so that par- 
ticles falling into that region are lost. This is a simplified 
but effective model for chaotic cavities with leads like the 
ones used in experiments with quantum dots. Quantum 
mechanically, the introduction of the hole corresponds to 
setting M rows (or columns) in the quantum propagator 
U to zero. Since it is no longer unitary, the new matrix 
U has left and right eigenstates 

C7|f^)=Z„|vI/«), (M/^|C/ = z„(vI/^| (1) 

and we may choose the following normalization and or- 
thogonality conditions: 

(*^1*^> = (*^l*^>, (*^|*^>=5„™. (2) 
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FIG. 1: Husimi representation of symmetrized right scar func- 
tions corresponding to a period 3 orbit of the triadic baker 
map, at iV = 81 (a) and = 243 (b). In panel (b) white 
crosses show the location of periodic points of the trajectory, 
and reflection by the diagonal produces a symmetric partner. 
Absolute value grows from white to black. 




eigenvalue number 

FIG. 2: In (a) we show the spectrum of the open triadic 
baker map for A'^ = 81 (circles) and the spectrum of the scar 
matrix (crosses). Their moduli (ordered by decreasing value) 
are displayed in panel (b). 



The eigenvalues lie in the unit disk, = e^^" < 1, 

and the quantityTn, > is interpreted as the decay rate. 
It was shown in Q that in the semiclassical hmit the long- 
hved left and right eigenstates localize on the stable and 
unstable manifolds, respectively, of the classical trapped 
set (strange repeller). 

We associate with every primitive periodic orbit 7 of 
period L a total of L scar functions, as was done in [15| 
for closed maps. One starts by building what is called a 
tube function, or periodic orbit mode. 
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Ne,)}\q„p,), (3) 
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which is a linear combination of coherent states on each 
of the L distinct points (qj^Pj) of the orbit. Here 
9j ~ J2i=o where Si is the action acquired by the Ith 
coherent state in one step of the map. The total action 
of the orbit is eL = S^. The quantity = {NS^ + k)/L 
is a Bohr-Sommerfeld-like eigenvalue. 
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We denote by A the Lyapounov exponent of the system. 

The right and left scar functions associated with the 
periodic orbit are defined through the propagation, under 



the open map, of the tubes until around the Ehrenfest 
time Te = j\nN. Namely, 



^;^[7*e"-^!;*cos 



and 



/)^|C/*e-2"^'*cos 



(5) 



(6) 



The constants Af^'^ are chosen such that {ipj^kl'^j^k) = 
{ip!^^k\'^!^^k) and (V'^^fclV'^^fe) = 1. The cosine is used to 
introduce a smooth cutoff, and the propagation time t 
is taken of the order of Te- In contrast to what is done 
for closed systems, we do not use negative powers of the 
matrix U, i.e. the tubes are propagated in only one direc- 
tion in time. This implies that the phase space support of 
right and left scar functions becomes localized on the un- 
stable and stable manifolds, respectively, of the periodic 
orbit, in consonance with the properties of resonances. It 
is natural to order these resonant scar functions accord- 
ing to the modulus of {fp^ kl^l'^jk) that longest-living 
ones come first. Finally, we note that it is convenient to 
impose the symmetries of the map on these functions (if 
two orbits are related by symmetry, we build symmetric 
and antisymmetric scar functions). 
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FIG. 3: Husimi representation of right resonances number 3 
and 12 (ordered by decreasing eigenvalue moduli) correspond- 
ing to the triadic baker map at A*' = 81. 



FIG. 4: Husimi representation of two symmetrized right scar 
functions corresponding to the same orbit of period 5, at = 
81. White crosses show the periodic points of the orbit, and 
reflection througli the diagonal produces a symmetric partner. 



It is by now established [IJ, [IJ, |15| that it is only nec- 
essary to use short periodic orbits to obtain good approx- 
imations to quantum spectra and eigenstates. By 'short' 
we mean orbits with periods up to around the Ehren- 
fest time of the system. This is not unexpected, because 
Te is the time when quantum interference effects become 
important. In the present case, since all periodic points 
are on the trapped set, the theory approximates only the 
long-lived states. Using the ordering mentioned above we 
construct the matrix (■!/j^|/7|'0^), which we call the scar 
matrix, as an approximation to the 'long-lived sector' of 
U. This matrix is by construction almost diagonal (in the 
sense that it is equal to a diagonal matrix plus a sparse 
one). 

What is the dimension of the scar matrix? For chaotic 
systems the number of periodic points grows with the pe- 
riod L like e'^^ where h is the topological entropy. Tak- 
ing orbits with periods up to Te we have e'*"^^ periodic 
points and corresponding scar functions. However, for 
small openings h is related to the fractal (information) 
dimension of the trapped set by d = 2h/X Since 
e^'^'^ = N we conclude that the matrix dimension (and 
the number of long-lived states) scales with N as N'^^^, in 
agreement with the fractal Weyl law 0| . Note that the 
dimension of C/ is TV, so our approach leads to a power 
saving in the size of the matrices used. This therefore 
represents a considerable improvement in efficiency. 

The above reasoning is in a sense complementary to 
the one presented in [<y\. There the authors considered 



quantum states which escape from the system before the 
Ehrenfest time (short-lived states). As a consequence 
they were led to regions of phase space that are preimages 
of the hole. Conversely, we are attempting to construct 
the quantum states which do not escape from the system 
before the Ehrenfest time (long-lived states) and are thus 
lead to short periodic orbits on the repeller. Consistently, 
both approaches result in the fractal Weyl law. 

As an example of the formalism, we use the triadic 
baker map, as in |8j. For this map the Lyapounov expo- 
nent is A = In 3, and we choose the quantum dimension 
to be TV = S'^ so that Te = k. The trapped set is the 
cartesian product Can x Can where Can is the usual 
middle-third Cantor set of dimension In 2/ In 3. The frac- 
tal Weyl law therefore predicts that the number of long- 
lived states should grow like A^'" 2/ in 3 _ 2^ (^^ig jg actu- 
ally the exact number if Walsh quantization Q is used). 
Let us take for instance fc = 4 and build scar functions for 
orbits with period up to 5 (there are 51 periodic points 
in total). We illustrate this construction in Figure [TJ 
where we show the Husimi plots of a symmetrized right 
scar function corresponding to an orbit of period 3, at 
= 81 and N — 243. It can clearly be seen that the 
probability extends along the unstable manifolds of this 
periodic orbit {q axis direction). We have used t = Te- 

In Figure [2] we present the exact quantum spectrum 
and the spectrum of the scar matrix (solution of a gen- 
eralized eigenvalue problem since {il'n\'4'm) 5nm)^ both 
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for TV = 81. We see that the latter provides excellent 
approximations to the first 30 resonances: they are all 
reproduced accurately and moreover there are no spuri- 
ous eigenvalues among them. We have verified that for 
A'^ = 3^ and using orbits up to period 6 (matrix dimension 
106) the first 55 resonances are reproduced accurately 
and without spurious eigenvalues. We have also verified 
that the method works well for other baker maps. 

Figure [3] shows the Husimi functions of the right reso- 
nances number 3 and 12 (ordered according to decreasing 
moduli of eigenvalues), again for = 81. They are both 
strongly localized around periodic orbits, i.e. they are 
scarred. The corresponding eigenstates of the scar ma- 
trix are indistinguishable from the exact ones, showing 
that this matrix is indeed a good approximation to the 
long-lived sector of U. 

In Figure H] we show two symmetrized scar functions 
built from the same periodic orbit, of period 5. Because 
some of the periodic points are very near the symme- 
try lines in phase space (the diagonals of the square) , we 
see that interference between the orbit and its symmetric 
partner makes the functions look rather different. Note 
the similarity between Figure a) and the scar function 
of Figure HJa) . A single element of our base captures 
almost all the structure of an exact quantum eigenstate. 
On the other hand, the resonance shown in Figure 
results essentially from the combination of the scar func- 
tion in Figure [ija) and the one in Figure lUJb) . 



We postpone a more detailed analysis to a future publi- 
cation, but the convenience of our approach to the study 
of the phase space morphology of quantum resonances is 
clear. Indeed, scar functions have proved extremely use- 
ful in the study of scarring effects, providing for example 
ways to quantify scarring and to understand the influence 
of homoclinic motion on scars fl^ . We expect it will also 
permit a more systematic study of scarring effects in open 
systems a subject still in its infancy. For instance, 
there is certainly an interesting interplay between scar- 
ring and the decay rate, so that more scarred states are 
expected to live longer. We believe our approach will 
shed some light on this issue. 

To conclude, we have introduced a theory based on 
short periodic orbits for quantum chaotic scattering. It 
may be argued that the use of the scar matrix offers no 
real advantage since its construction makes use of the 
quantum propagator. However, previous studies of closed 



systems [IJ] suggest that semiclassical approaches can 



be successfully implemented within this framework, be- 
cause the propagation times involved are not longer than 
Te- We are currently working in this direction. An- 
other direction to follow is to adapt the theory to dielec- 
tric boundary conditions in order to treat microlasers, an 
application where spectacular manifestations of scarring 
can be observed 0]. 
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